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The first moon tracking computer program included here was 
rewritten from WA1JXN/WA3GPL and K5JL versions presented in earlier 
issues of the EME notes. Warren Butler, W2WD, wrote the program for 
TRS-80 level II BAS I C ( 16K) . Cassette copies of the program are 
available from Warren for the cost of the cassette and postage 
(approximately $1.00). 

The second program was written in FORTRAN IV by Geoffrey Grayer, 
G3NAQ. Again the WA1JXN/WA3GPL program served as the starting point 
for this effort. 
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COORDINATES OF THE MOON 


10 1 

20 PRINT 

30 * DISLAYS GHA, DECLINATION, AZIMUTH, ELEVATION OF THE MOON 
40 * FROM A SELECTED LATITUDE, LONGITUDE FOR SELECTED DATES 
50 • AND TIMES (GMT). EASTERN STANDARD TIME IS SHOWN BUT CAN BE CHANGED— 
60 * @ LINE 2040, CHANGE 500 TO TIME DIFFERENTIAL DESIRED— 

70 1 EG., FOR PACIFIC STANDARD TIME, USE 800 IN PLACE OF 500. 

80 * @ LINES 850 S 860 CHANGE EST TO TIME ZONE DESIRED— 

90 • EG., REPLACE EST WITH PST FOR THE EXAMPLE BEING USED. 

100 1 HARDCOPY OUTPUT CAN BE SELECTED IF PRINTER IS AVAILABLE; 

110 * OTHERWISE, DATA WILL BE DISPLAYED ON CRT ONLY. 

120 PRINT 

130 1 BASED ON PROGRAMS BY LANCE COLLISTER, WA1 JXN/WA3GPL AND 
140 * JAY LI E EM ANN, K5JL. 

150 ' MODIFIED FOR TRS-80 LEVEL II BASIC (16K) BY WARREN BUTLER, W2WD 
160 PRINT (4/4/79). 

170 ' DATA FOR UP TO 31 DAYS CAM BE REQUESTED EACH TIME THE 

180 1 COMPUTER ASKS FOR AN INPUT BY PRINTING ?, ENTER DATA IN 
190 ' THE FORMAT REQUESTED. AFTER THE LAST INPUT, INSERT ZEROS 
200 * IN THE SAME FORMAT. 

210 1 ANSWER OTHER QUESTIONS AS APPROPRIATE 

220 * UNIVERSAL WINDOWS ARE SHOWN BY LETTERS FOLLOWING DEC VALUES 

230 1 U = EUROPEAN UNIVERSAL WINDOW 

240 1 W = W/VE UNIVERSAL WINDOW 

250 ' J = J/VK/ZL UNIVERSAL WINDOW 

260 CLEAR 500 

270 DIM F(31),V(31),Y(31),Q(31),S(31) 

280 P5=2. 0000000000-3. 1415926535 

290 D5 =360. 0000 0000000/P5 
300 R5=P5/360. 0000000000 

310 CLS : PRINT : PRINT : PRINT :PRINT"WHAT ARE THE STATION CALL LETTERS"; 

320 INPUT W$ 

330 PR I NT"WHAT IS YOUR LATITUDE IN DEGREES, MINUTES"; 

340 INPUT L5,U5 

350 PRINT"WHAT IS YOUR LONGITUDE IN DEGREES, MINUTES"; 

360 INPUT L6,U6 
370 L5=(L5+U5/60)»R5 

380 L6 = (L6+U6/60) :: R5 

390 PRINT"WHAT IS THE DESIRED PRINTING INCREMENT IN MINUTES"; 

400 INPUT I 

410 I NPUT"DO YOU ONLY WANT PRINTOUT WHEN THE MOON IS NEAR THE HORIZON 

0)";B$ (YES/ 

430 IF B$="YES" THEN 460 
440 LET 16=100 
450 GOTO 550 

460 I NPUT"BE LOW WHAT ELEVATION IN DEGREES DO YOU WANT PRINTOUT"; 16 
470 I NPUT"DO YOU WANT HARDCOPY PRINTOUT (YES/NO)"; WW$ 

480 PRINT"WHAT ARE THE GMT MONTH, DAY, YEAR DESIRED?" 

490 PR INT"USE FORMAT Mi1, DD, YYYY 4-DIGITS FOR YEAR" 

500 FOR N=1 TO 31 
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51 0 INPUT F(lO, VGO, Y GO 
520 IF F 00=0 THEN 640 
530 NEXT N 
540 GOTO 500 

550 INPUT"DO YOU WANT HARDCOPY C YES/NO)" ; WW$ 

560 PRINT 

570 PR I NT"WHAT ARE THE GMT MONTH, DAY, YEAR, TIME BEGINNING, TIME ENDING?" 

580 PR INT"USE THE FORMAT MM, DD, YYYY, TTTT, TTTT 

REMINDER USE 4 DIGITS FOR YEAR!" 

590 FOR N=1T031 

600 INPUT fcn), voo, y(n), q(n), SOO 

610 IF F (N)=0 THEN 640 

620 NEXT N 

630 GOTO 590 

640 N5=N- 1 

650 FOR N=1 TO N5 

660 IF B$="YES" THEN 680 

670 GOTO 710 

680 E 1 = 2400 

690 B=0 

700 GOTO 730 

710 E1=S(N) 

720 B=Q(N) 

730 M=FCN) 

740 D=V(N) 

750 Y=Y(N) 

760 Y 1= Y-( INTCY / 100) **100) 

770 PRINT 

780 IF WW$ = "YE S" LPRINT 

790 PR I NT : C LS 

800 IF WW$="YES" LPRINT 

810 PRINT "POSITI ON OF THE MOON ON " ;M; "/" ; D ; " ; Y 1 ; " GMT FROM"" "W$ 

820 IF WW$=" YES" LPRINT "POSITI ON OF THE MOON ON ";M; "/" ; D; "/" ; Y 1 ; " GMT FROM"" "W 

$ 

830 PRINT 

840 IF WW$ = " YE S" LPRINT 

850 PRINT"GMT"TABC11) m GHA"TABC21) ,, DEC"TABC35)"EST"TA3C47)"AZ"TABC57) m EL" 

860 IF WW$ =" YES" LPR I NT"GMT" TABC 1 1 )"GHA"TAB( 2 l ) " DEC"TAB( 35 )"ES T" TAB$ 47 )"AZ" TAo( 5 
7)"EL" 

870 PRINT 

880 IF WW$ =" YE S" LPRII>JT 
890 11=2 

900 IF M>=3 THEN 980 

910 IF INT(CY-1853)/4)<11 THEN 940 

920 Cl = -1 

930 GOTO 950 

940 C 1=0 

950 Jl = 365 :: CY-185 3) + D+30 :; CM+9)+INTCCM+10)/2) 

960 J2=INTCCY-1853)/4)+l+Cl 
970 GOTO 1090 

980 IF INTCCY-1852)/4XH THEN 1010 
990 Cl=-1 
1000 GOTO 1020 
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10 10 C1=0 

1020 IF M=9 THEN 1060 
1030 IF M-ll THEN 1060 
1040 C2-0 
1050 GOTO 1070 
1060 C2 = 1 

1070 Jl = 365 :: CY-185 2) + D+30 ;: CM-3)+lNTCCM-2)/2) 

1080 J2=INT(CY-1852)/4)+Cl+C2 

1090 J=J 1+ J2 
1100 Tl= J-l 7472.5 

1110 D9=CB-INTCB/100)-100)+INTCB/100)»60 

1120 D6=CE1-INT(E1/100)«100)+INT(E1/100>“60 

1130 D7=D9-D6 

1140 D8=D7-I 

1150 IF D 7>0 THEN 1170 

1160 GOTO 1190 

1170 IF D8> =0 THEN 2150 

1180 B = E1 

1190 T=(B-INT(B/100)::i00)/1440+INT(B/100)/24 
1200 T5=T1+T 

1210 K1 = (C . 75 12 13+.036601 102 :: T5)-INT(.751213 + .03660 U02«T5))i:p5 
12 20 K2 = ((,822513+. 0362916457 :: T5)-INT(.322513+. 0362916457 ;{ T5}) :: P5 

12 30 K3 = CC .995 766 + . 0027377785 2 :: T5)-INTC. 995 766 + .00273777852 :: T5)) :; P5 
12 40 K4 = CC.974271+. 0338631922 :: T5)-INTC. 974271+. 0338631922 “T5)) :: P5 
12 50 K5 = CC.0312525+. 036748195 7 :: T5)-I NT C.0312525 + . 0367481957 :: T 5)) :: P5 
12 60 L8 = Kl + .6 5 8”R5 :: SINC2”K4) + 6. 2 89 :: R5 :: SINCK2) 

12 70 L8 = L3-1.274 :: R5 :: SIN(K2-2 :: K4)-. 186 :: R5 ;c SIN(K3) 

12 80 L8 = L8 + . 2 1 4 :: R5 :: S I N (2 :: K2 ) - . 114 ;: R5 :: SIN(2 :: K5) 

12 90 L8 = L8-.05 9 :: R5 :: SIN(2 :: K2-2 :: K4)-.05 7 :: R5 :c SIN(K2+K3-2»K4) 

1300 K6 = K5 + .6593 :: R5 :: SIN(2 ;: K4) + 6. 2303 :; R5 :: SIN(K2)-1. 272 1: R5 :: SIN(K2-2”K4) 
1310 L7 = 5. 144”R5 :: SINCK6)-. 146 : 'R5 :: SINCK5-2«K4) 

1320 LET Dl = COS(L79 :: SIN(L8) :: .397821 + SlN(L7) :: .91746 3 
1330 LET D1=ATNCD1/CSQRC1-D1±2))) 

1340 Gl = 50+.5+CCDl)/C .792 )) :: D5 
1350 G2=80+CCD1)/(.808))»D5 
1360 G3=141.5-CCD1)»C.738)«D5) 

13 70 G4=170.5-CCD1)»C. 85 7) :: D5) 

1380 A2 = COS(L7) :: COSCL8)/COS(Dl) 

1390 Al = (COS(L7) :: SIN(L8) :: . 9 17463-S I N(L7) . 39 78 2 1 ) /COS(Dl ) 

1400 A=ATN(A1/A2) 

1410 GOSUB 1670 
1420 R 1=A 

1430 Ll=. 065709822 :: T1 

1440 L=T«2 4»1.002 7 38-f6.64605 5 + CLl-INTCLl/24)»24) 

1450 L=CL- INTCL/24) 5: 24) 

1460 G=(L/24) :: P5-R1 
1470 IF G4P5 THEN 1500 
1480 G=G-P5 
1490 GOTO 1530 
1500 IF G<0 THEN 1520 
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1510 GOTO 1530 
1520 G =G + P5 
1530 H=L6-G 

1540 E3 = COSCL5) :: COS(H)”COSCD1) + SINCdD : '*SINCL 5) 
1550 E2 = SQRCl-CE3 :c E3)) 

1560 E=ATN(CE3/E2)-C1/C61.3 3 :c E2))) 

1570 F=ATN(E3/E2) 

1580 IF E<0 THEN 2100 

1590 IF E*I6»R5 THEN 2100 

1600 A2 = sinCdO/(cosCl5) :: cosCf)) 

1610 A2=A2-CSirKL5)/COSCL5)) :: (SINCF)/COSCF)) 

16 20 A1 = SINCL5) :; SINCdO + COSCL5) :: COSCD1)”COSCH) 
1630 Al=CSINCH)»COS(Dl))/SQRCl-Al±2) 

1640 A=ATN(A1/A2) 

1650 GOSUB 1670 

1660 GOTO 1820 

1670 IF A=0 THEN 1690 

1680 GOTO 1730 

1690 IF A2<0 THEN 1710 

1700 GOTO 1810 

1710 A=P5/2 

1720 GOTO 1810 

1730 IF A>0 THEN 1790 

1740 IF A2<0 THEN 1770 

1750 A=P5+A 

1760 GOTO 1810 

1770 A=P5+(A-P5/2) 

1780 GOTO 1810 

1790 IF A2=>0 THEN 1810 

1800 A= A+P 5/2 

1810 RETURN 

1820 IF (T-I 1)>C2 :: I)/1440 THEN 1840 
1830 GOTO 1850 
1840 PRINT 

1850 BS=INT(3+.5):BS$= ,, ~ :: H" 

1860 Zl = INT(A :: D5 :c 10+.5)/10 
1870 Z2 = INTCE»D5 :: 10 + . 5)/10 
1880 Z3=INT(G :: D5 :: 10+.5)/10 
1890 Z4=INT(Dl :: D5 :: 10+.5)/10 
1900 IF Z4<0 THEN 2030 
1910 IF Z3<G1 THEN 2030 

1920 IF Z3>G2 THEN 1940 

1930 GOTO 1970 
1940 IF Z3<G3 THEN 1990 

1950 IF Z 3>G4 THEN 2030 

1960 GOTO 2010 
1970 Y$="U" 

1980 GOTO 2040 
1990 Y$="W" 

2000 GOTO 2040 
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20 10 Y $ = *' J" 

2020 GOTO 2040 
2030 Y$ = " f ' 

2040 ES=(INT(B+.5))-500 
2030 IFES<=0THENES=ES+2400 
2060 ES$ = " — H" 

2070 PRINTUS INGBS$; BS; : PR I NTTABC 1 0)Z3TAB(20) Z4; Y $TAB( 33 )US I NGES$ ; ES; : PR INTTAB$ 4 
) Z 1TABC 5 5 ) Z2 

2080 IF WW$ = ,, YES M LPR I NTUS I NGBS $ ; BS ; *. LPR I NTTAB( 10) Z3TAB( 20 ) Z4; Y$TAB <3 5 )US I NGE S $ ; 

ES; :LPRINTTAB(45)Z1TAB(55)Z2 
2090 I 1=T 
2100 B=B+ I 

2110 Z=(B-INT(B/100)”100)-60 

2120 IF Z<0 THEM 1110 

2130 B=INTCB/100)-100+100+Z 

2 0 GOTO 1110 

2150 NEXT N 

2160 N=0 

2170 PRINT 

2180 PRINT 

2190 PRINT"DO YOU WANT MORE INFORMATION (YES/NO)"; 

2200 INPUT D$ 

2210 IF D$="YES" THEN 280 
2 ° °0 END 
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POSITION OF THE MOON ON 5 / 19/79 GMT FROM W 2 WD 


GMT 

GHA 

DEC EDT 

AZ 

EL 

”615 

357.4 

-10 

215 

106 

2.2 


1 

- 9.9 

230 

108.4 

4.9 

”645 

4.6 

- 9.9 

245 

110.8 

7.5 


8.2 

- 9.8 

300 

113.4 

10.1 

"715 

11.8 

- 9.8 

315 

115.9 

12.6 

”730 

15.4 

- 9.8 

330 

118.6 

15.1 

”745 

19.1 

- 9.7 

345 

121.3 

17.5 

”800 

22.7 

- 9.7 

400 

124.2 

19.8 

”815 

26.3 

- 9.6 

415 

127.2 

22.1 

”830 

29.9 

- 9.6 

430 

130.3 

24.3 

”845 

33.5 

- 9.6 

445 

133.5 

26.4 

”900 

37.1 

- 9.5 

500 

136.9 

28.4 

”915 

40.8 

- 9.5 

515 

140.4 

30.3 

”9 30 

44.4 

- 9.4 

530 

144.1 

32 

”945 

48 

- 9.4 

545 

147.9 

33.6 

1000 

51.6 

- 9.4 

600 

151.9 

35 

1015 

55.2 

- 9.3 

615 

156.1 

36.2 

1030 

58.8 

- 9.3 

630 

160.4 

37.3 

1045 

62.5 

- 9.2 

645 

164.8 

38.2 

1100 

66.1 

- 9.2 

700 

169.3 

38.8 

1115 

69.7 

- 9.1 

715 

174 

39.3 

1130 

73.3 

- 9.1 

730 

178.6 

39.5 

1145 

76.9 

- 9.1 

745 

183.3 

39.5 

1200 

80.5 

-9 

800 

188 

39.3 

1215 

84.2 

-9 

815 

192.6 

38.8 

1230 

87.8 

- 8.9 

830 

197.1 

38.2 

1245 

91.4 

- 8.9 

845 

201.5 

37.3 

1300 

95 

- 8.9 

900 

205.8 

36.2 

1315 

98.6 

- 8.8 

915 

210 

35 

1330 

102.2 

- 8.8 

930 

214 

33.5 

1345 

105.9 

- 8.7 

945 

217.9 

31.9 

1400 

109.5 

- 8.7 1000 

221.5 

30.2 

1415 

113.1 

- 8.6 1015 

225.1 

28.3 

1430 

116.7 

- 8.6 1030 

228.5 

26.4 

1445 

120.3 

- 8.6 1045 

231.7 

24.3 

1500 

123.9 

- 8.5 1100 

234.8 

22.1 

1515 

127.6 

- 8.5 1115 

237.8 

19.8 

1530 

131.2 

- 8.4 1130 

240.7 

17.5 

1545 

134.8 

- 8.4 1145 

243.5 

15.1 

1600 

138.4 

- 8.4 1200 

246.2 

12.6 

1615 

142 

- 8.3 1215 

248.8 

10.1 

1630 

145.6 

- 8.3 1230 

251.4 

7.6 

1645 

149.3 

- 8.2 1245 

253.9 

4.9 

1700 

152.9 

- 8.2 1300 

256.3 

2.3 
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FORTRAN lv G I RELEASE 2.0 


MAIN 


DATE * 78111 


17/36/51 


PAGE 0001 


U0U1 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 

0012 

0013 

0014 

0015 

0016 


001 / 
001 8 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 


THIS PROGRAM IS DESIGNED TO CALCULATE THE AZIMUTH AND ELEVATION OF THE MOON 

THIS IS A FORTRAN VERSION OF PROGRAM 3 OF THE EIMAC EME NOTE AS-49-6 
ORIGINALLY WRITTEN BY LANCE COLHSTER WA3GPL/ CONVERTED BY GEOFF GRAYER G3NAQ. 

THE INPUT FORMAT IS AS FOLLOWS... 

FIRST DATA CARD: 

OPTIONAL IDENTIFICATION (AlU)/ LAT. OF STATION BEGS. (110)/ LAT. OF STATION 
MINS. (110)/ LONG. OF STATION DEGS. (110)/ LONG. OF STATION MINS. (110). 

SECOND AND SUBSEQUENT DATA CARDS: 

YEAR (110)/ MONTH (110)/ DAY (110)/ START HR S . MI NS . ( 1 1 0 ) / STOP H R S . M I NS . ( 1 1 0 ) 
I1ME INCREMENT MINS. (110) - DEFAULT 10 MINS./ M A X 1 MOM ELEVATION DEGS. (110) - 
MAY BE OSED TO SELECT PRINTOUT ONLY WHEN THE MOON IS NEAR THE HORIZON - 
DEFAOLT 90 DEGS. 

ALL TIMES ARE IN GMT OSING THE UOUO TO 2400 HOUR SYSTEM. 

A BLANK CARD DENOTES THE END OF THE D A 1 A SET AND TERMINATES THE PROGRAM. 
PRINTOUT IS SUSPENDED WHENEVER ThE ELEVATION OF THE MOON IS NEGATIVE. 

FOR FURTHER INFORMATION/ REFER TO THE EIMAC NOTE. 

GEOFFREY H. GRAYER G3NAQ 
BRIGHTWALTON/ BERKSHIRE 
APRIL/ 1978 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

INTEGER D/Y/C1/C2/B/Z/E1/DO/D7/D8/DV/GMT 

f* I/P5/D5/R5/K1 / K2/K3/K4/K5/L8/K6/L7/D1 / A1 / A2/A/L1 / L/ E3/E2 
KbAL Lb#L6#I1 
R E AL*4 JU 

DIMENSION IDENT (1U)/GMT(4) 

FNA(X)=AlNr(X*D5*1U.+U.5)/10. 

FNB(X)=(X-AINT (X) ) *P 5 
PI=3.14l5v26535 
P5*2 . *P I 
05=180. /PI 
R5 = PI/iai). 

C READ STATION IDENT/ LAT DEGS/ LAT Ml NS/ LONG DEGS/ LONG MINS 
READ (5/1) I DEnT /LATD/LA TM/LONgD/LONGM 

1 FORMAT (1UA1/4I1U) 

L5=(LATD + LATM/60. ) *R5 
L6=(L0NGD+L0NGM/6U.)*H5 

10 CONTINUE 

C ^^ A (MAXlMUM°fcLEVATION S |)EGS bl1,, ' (feND T IMfc G*T),(TIME INCREMENT MINS)/ 

READ (5/2) Y/M/D/d/EI/I/16 

2 FORMAT (7110) 

IF(Y.EQ.U) STOP 
IF(I.EQ.U) 1*10 

1 F ( I6.EU.U) 16 = 100 
IF(EI.Eu.O) E 1 =24 Ou 
11=2 

C CALCULATE JULIAN DATE 

IF (M.GE.3) GO TO 16 
C1=-1 

1F(((Y“1853)/4).LT.11) C1=0 

J 1 =365* ( r-1 853 ) ♦ 3U*(M+9) ♦ ( (M + 1 0) /2) ♦ D 
J2=( ( Y-1853) /4) + Cl ♦ 1 
GO TO 2/ 

16 C1=-1 

IF (( (Y-1 852 )/4).LT.11) C1=U 
IF (M.E0.9.0R.M.EQ.1 1 ) Gu TO 24 
22 C 2 = 0 

GO TO 25 

24 C 2 = 1 

25 J 1 = 365* ( Y-1 852 ) ♦ 3u*(M-3) + ( ( M— 2 ) / 2 ) + D 
J 2= ((Y-1852)/4) + Cl ♦ C2 

27 J * J 1 -f j 2 
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FORTRAN IV G1 RELEASE 2.0 


MAIN 


DATE 


78111 


17/56/51 


PAGE 0002 


0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 

U052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 
006 7 
0068 

0069 

0070 

0071 

0072 

0073 

00/4 

0075 

0076 

0077 

0078 

0079 

0080 

0081 

0082 

0083 

0084 

0085 

0086 

008/ 

0088 

0089 

0090 

0091 

0092 

0093 

0094 

0095 

0096 
009/ 

0098 

0099 

0100 
0101 
0102 

0103 

0104 

0105 


LATS '#•214, 
1 * FI 2 .1 / / ) 


T1®FL0AT(J)-17472.5 
J D® FLOAT (J)+2397547.5 

WRITE ( 6# 3) IDENT,LATD,LATM,LONGD,LONGM,Y,M,D, J D 
3 FORMAT ( * 1 LON AR COORDINATES FOR ',10A1/' STATION 
♦ LONG: ' ,214/ ' DATE ' , 14/ ' / ' , 1 2, ' / 1 / 1 2 , ' JD< 

WRI TE (6,5 ) 

5 FORMAT ( ' GMT AZ EL GHA DEC'/) 

29 D9*B - ( B / 1 00 ) *1 00 ♦ (B/1UU)*6U 

06® El - ( El / 1 00 ) *1 00 ♦ <E1/1U0)*60 

D7®D9-D6 

D8»D7-I 

IF (D7.LE.0) GO TO 38 
IF (D8.GE.0) GO TO 10 
B*E 1 

C CALCOLATE LONAR LAT AND LONG 

38 T*FLOAT(8-(B/10U)*100)/144U.+FLOAT(B/100) / 24. 

T5«T1+T 

K1 ®FNB ( 0.Z51 21 3+0. 0366011 02*T 5) 

K2-FNB (0.82251 3+0. 036291 64 5 /*T5) 

K3*F NBC 0.995766 + 0. 002 7377 7852* T5) 

K 4® F NBC 0.9742 /I +0.0338631 92 2*T5) 

K 5® F NB ( 0.031 2 5 25+0. 036/481 95 /*T 5) 
L8®Kl+0.6S8*R5*SlN(2.*»U)+6.289*R5*SIN(K2) 

L8®L8-1 -2 74*R5*S1N(K2-2.*K.4)-U.186*R5*SIN(K3) 

L8»L8+0.214*R5*SIN(2.*K2)-0.114*R5*SINC2.*K5) 

L8®L8-0.059*R5*S1N(2.*K2+2.*K4)-U.057*R5*SIN(K2+K3-2.*K.4> 

K6*K5+.6593*R5*SIN(2,*K4)+6.23U3*R5*SIN(k2)-1.272*R5*SIN(k2-2.*K4) 

L7«5.144*R5*SIN(K6)-0.146*R5*SIN(K5-2.*K4) 

C CALCOLATION OF RA AND DEC „ „ . 

D1®C0S(L 7 )*SIN(L8)*0.39/821+SIN(L/)*0. 917463 
D1®ATAN2(D1 # SORT (1 .-D1 **2) ) 

A2*CQS(L/)*COS(L8)/COS(Dl) 

A1 « (COS (L/)*S IN (L8) *0.91 7463-SIN (L 7) *0.39/821 )/ COS (D1 ) 

A«A T AN2 ( A1 , A2 ) 

R 1 * A 

L1®0. 065/09822*11 

L*T *24 . *1 .00 2/38 + 6. 646055+(Ll-AlNT (Li /24)*24.) 

L*(L-AINT (L/24)*24) 

C CALCOLATION OF GREENWICH HOUR ANGLE G FROM LOCAL S 1 DERI AL TIME 
G* ( L / 24 ) *P5-R 1 
IF(G.LT.P5) GO TO 67 
G*G-P 5 
GO TO 71 

67 IF(G.LT.O.) GO TO 69 
GO TO 71 
69 G**G + P5 

C CALCULATION OF LOCAL HOUR ANGLE H FROM GHA 
71 H»L6-G 

C CALCULATION OF ELEVATION E 

E3«C0S(L5)*C0S(H)*C0S(D1)+SIN(D1)*SIN(L5) 

E2«SQRT(1-E3**2) 

£«AT AN2 ( E3,E2 ) 

IF (E.LT.O.) GO TO 117 
IF (E.GT. ( 1 6* R 5 ) ) GO TO 11/ 

C CALCULATION OF AZIMUTH A 

A2»SIN(D1 ) / (COS (L5)*C OSCE) ) 

A2*A2-(SIN(L5)/C0S(L5))*(S1N(E)/C0S(E)) 

A1®SIN(L5)*SIN(D1)+COS(L5)*COS(D1)*COS(H) 

A1*(SIN(H)*C0S(Dl ) ) / SORT (1 .-A1**2) 

A-ATAN2CA1 ,A2) 

AZ*FNA( A) 

EL=FNA(E) 

GH A* F N A ( G ) 

DEC»FNA ( D1 ) 

GMT ( 1 ) *B / 1 UOU 
GMT (2)®B/1U0— (B/10U0)*1U 
GMT (3)»B/1U-(B/100)*1 0 
&MT(4)*B-(B/1U)*1U 

IF ( (T-I1 ) .GT. (2. *1/1440.) > WRITE (6,7) 

/ FORMAT (1H ) 

1U4 WRITE(6,6) GMT, AZ, EL, GHA, DEC 

6 FORMAT (3H ,4I1,4F12.1) 

1 1 ■ T 

11/ B®B+I 


FORTRAN IV G1 RELEASE 2.0 


0106 

0107 

8181 

0110 


MAIN 


Z®B°(B/100)*10Q - 60 
IF(Z.LT.O) GO TO 29 
8®<B/100)*100 ♦ 100 ♦ Z 
GO TO 29 
END 


DATE ■ 78111 17/56/51 PAGE 0003 
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LUNAR COORDINATES FOR G3NAQ 


AT ION LAT; 51 34 

LONG: 

1 18 


ITE 197 a/ 

5/20 J D* 

2443648.0 



GMT 

A2- 

EL 

GHA 

DEC 

0000 

-143.7 

24.3 

34.2 

-7.7 

0010 

-141.3 

23.3 

36.6 

-7.7 

0020 

-158.9 

22.4 

39. U 

-7.7 

0030 

-136.6 

21.3 

41.5 

-7.8 

0040 

-154.4 

20.3 

43.8 

-7.8 

0050 

-152.1 

19.1 

46.3 

-7.8 

0100 

-129.9 

18.0 

48.7 

-7.8 

0110 

-127.7 

16.8 

51.1 

-7.9 

01 20 

-125.6 

15.6 

53.4 

-7.9 

0130 

-125.4 

14.3 

56.0 

-7.9 

0140 

-121.4 

13.0 

58.4 

-7.9 

0150 

-119.4 

11.7 

60.7 

-8.0 

0200 

-117.4 

10.4 

63.1 

-8.U 

0210 

-113.4 

9.0 

65.5 

-8.0 

0220 

-113.3 

7.6 

68.0 

-8.0 

0230 

-111.5 

6.2 

70.4 

-8.1 

U240 

-1 09.6 

4.8 

72.8 

-8.1 

0250 

-107.6 

3.3 

75.2 

-8.1 

0300 

-105.8 

1.9 

77.6 

-8.2 

0310 

-103.9 

0.4 

80.0 

-8.2 

1 720 

107.2 

0.1 

285.0 

-10.5 

1 7iU 

109.1 

1.5 

287.4 

-10.5 

1740 

111.0 

2.8 

289.7 

- 1 0.5 

1 750 

1 15.0 

4.2 

292.1 

-10.6 

1800 

114.9 

5.5 

294.6 

-10.6 

1810 

1 16.9 

6.9 

29 7.1 

— '1 U • 6 

1 82u 

1 18.9 

8.1 

299.4 

-1 0.6 

185U 

120.9 

9.4 

301.8 

-10.7 

1840 

122.9 

10.7 

304.2 

-1U. 7 

1850 

123.0 

11 .9 

306.6 

-10.7 

1900 

127.1 

13.1 

309. U 

-1 o . 7 

1910 

129.2 

14.2 

311.4 

-10.8 

1920 

151 .4 

13.4 

313.9 

-10.8 

1950 

155.6 

16.5 

316.3 

- 1 u.a 

1 94U 

153.8 

1 7.5 

318.7 

-10.8 

1950 

158.1 

18.3 

321 .1 

-1 U.9 

2000 

140.4 

19.5 

323.5 

-10.9 

2010 

U2.7 

20.4 

325.9 

-10.9 

2020 

143.0 

21 .2 

328.3 

-10.9 

2030 

147.4 

22. u 

330.8 

-1 1 .U 

2040 

149.9 

22.8 

333.2 

-1 1 .0 

2050 

132.5 

23.5 

335.6 

-11 .0 

2100 

154.7 

24.1 

S3 7 . 9 

-11.1 

21 1 u 

137.3 

24.7 

34 0.4 

-1 1.1 

2120 

1 5 v . 9 

25.2 

342.8 

-11.1 

21 30 

162.4 

25.7 

345.1 

-11.1 

2140 

16 5. 0 

26.1 

34 7.5 

-1 1 .1 

21 50 

167. 7 

26.4 

350.1 

-11.1 

2200 

1 7U.3 

26.7 

352.5 

-11.2 

2210 

172.8 

26.9 

354.8 

-1 1 .2 

2220 

17 5.5 

27.0 

33 7.2 

-1 1.2 

2230 

178.2 

27.1 

359.7 

-11.2 

2240 

-1 79.1 

27.0 

2.0 

-1 1.3 

2250 

-176.5 

27.0 

4.4 

-11.3 

2300 

-178.8 

26.8 

6.8 

-11.3 

2310 

-171.2 

26.6 

9.3 

-1 1 .4 

2320 

- 1 68.6 

26.3 

11.7 

-11.4 

2330 

-166.0 

26.0 

14.1 

-11.4 

2340 

-163.4 

25.5 

16.5 

-11.4 

2550 

-160.6 

25.1 

18.9 

-1 1.5 

2400 

- 1 38.5 

24.5 

21.3 

-11.3 
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